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The recent proliferation of correlated percolation models — models where the addition of 
edges/vertices is no longer independent of other edges/ vertices — has been motivated by the quest 
to find discontinuous percolation transitions. The leader in this proliferation is what is known as 
explosive percolation. A recent proof demonstrates that a large class of explosive percolation-type 
models does not, in fact, exhibit a discontinuous transition[0. Riordan and L. Warnke, Science, 
333, 322 (2011)]. We, on the other hand, discuss several correlated percolation models, the fc-core 
model on random graphs, and the spiral and counter-balance models in two-dimensions, all exhibit- 
ing discontinuous transitions in an effort to identify the needed ingredients for such a transition. 
We then construct mixtures of these models to interpolate between a continuous transition and a 
(— i ' discontinuous transition to search for a tricritical point. Using a powerful rate equation approach, 

we demonstrate that a mixture of k — 2-core and k = 3-core vertices on the random graph exhibits a 
tricritical point. However, for a mixture of fc-core and counter-balance vertices, heuristic arguments 
and numerics suggest that there is a line of continuous transitions as the fraction of counter-balance 
vertices is increased from zero with the line ending at a discontinuous transition only when all 
vertices are counter-balance. Our results may have potential implications for glassy systems and a 
recent experiment on shearing a system of frictional particles to induce what is known as jamming. 

I. INTRODUCTION 

Percolation is the study of connected structures in disordered networks [lj. For example, two edges meeting at a 
vertex form a connected structure called a cluster of size two. As more edges are randomly and independently added 
to the network, the average cluster size grows until there ultimately exists a spanning cluster in finite-dimensional 
lattices (or a giant component in random graphs). The onset of a spanning cluster, which is indicative of a transition 
from a non-spanning to spanning phase, exhibits properties of a continuous phase transition. The simplicity of this 
nontrivial model allows one to catalog many of its properties such that it is the Ising model of geometrically-driven 
phase transitions @, Q . 

While the simplicity of percolation is part of its power, there has been a recent renaissance in developing models 
beyond ordinary percolation in an effort to discover new types of transitions such as a discontinuous one. The 
main driving force behind this endeavor is what is known as explosive percolation 4-8]. To be specific, two edges 
are considered at random and the edge that minimizes, for example, the product of the two clusters it joins is then 
retained and the other discarded, i.e. a choice has been made as to which edge to retain. Initial numerical data for this 
particular model on random graphs suggested that the percolation transition is discontinuous such that the emergence 
of the giant component is explosive, hence the term explosive percolation. Because explosive percolation goes beyond 
ordinary percolation where there is no choice between edges, a different type of transition is not necessarily surprising. 
. However, a recent mathematical proof shows that the transition in this class of models involving choice on random 
t— I ' graphs is, in fact, continuous [9l. flo| . 

Other explosive percolation-type models are being investigated for the possibility of a discontinuous transition. 
. ' For instance, simulations of explosive percolation on finite-dimensional lattices show signs of a discontinuous phase 
transition [ll|, [l2| . In addition, in light of the recent proof that models involving a choice (Achlitopas processes) 
H 1 on random graphs exhibit a continuous transition, some researchers have begun to study other models beyond or- 
dinary percolation where there exist various constraints on the occupation of edges (and/or vertices). For example, 
the Bohman-Frieze-Wormald model allows for the addition of an edge provided it participates in a cluster smaller 
than some prescribed size with the prescribed size being updated as edges are added [HI . Recent work suggests a 
discontinuous transition on random graphs for this model [lj, [l5[ and even more recent work suggests a similar result 
on finite-dimensional graphs such as the cubic lattice [l6j . 

While work progresses on these more complicated percolation models, perhaps the simplest beyond ordinary per- 
colation model is known as fc-core percolation [l7l - [T9j . fc-core percolation is defined as the following: Each vertex on 
a graph needs at least k occupied edges; if the constraint is not obeyed, the vertex and the edges attached to it are 
recursively removed until a stable fc-core configuration is reached (with every vertex obeying the fc-core constraint). 
It turns out that fc-core percolation on random graphs exhibits a continuous transition for fc < 2 and a discontinuous 
transition for k > 3. However, unlike a typical discontinuous transition there exists not one, but several, diverging 
lengthscales such that the transition is an unusual one. 

So there indeed exists discontinuous percolation transitions in mean-field. What about finite-dimensions? It turns 
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out that fc-core percolation on various finite-dimensional lattices, such as the triangular lattice [20|, either falls into 
the same universality class as ordinary percolation or does not exhibit a transition (but interesting finite-size effects) , 
i.e. p c = 1 [HlHSl, where p c denotes the critical occupation probability. Of course, in high enough dimensions it is 
conjectured that the fc > 3-core transition becomes discontinuous .23]. While obtaining a discontinuous transition 
with fc-core in finite-dimensions has been difficult, a new class of correlated percolation models, dubbed jamming 
percolation, with constraints more complex than fc-core, has been recently shown to exhibit a discontinuous transition 
in finite-dimensions (2~3 - [27j . In addition to exhibiting a discontinuous transition, these models also exhibit diverging 
lengthscales that grow either as a power law or faster than power law such that they, too, are not the garden variety 
discontinuous transition. 

Since there exists percolation models beyond ordinary percolation exhibiting a discontinuous transition, it is possible 
to construct hybrid models where some fraction of the edges are occupied via ordinary percolation and the remaining 
edges are occupied via a choice method or some other constrained method. Since ordinary percolation dominates in 
one limit and the correlated percolation model, which includes explosive percolation, dominates in the other limit, 
it may be possible to locate a tricritical point bordering the continuous and discontinuous regimes. The proof that 
explosive percolation is continuous on random graphs removes the possibility of a tricritical point on random graphs 
using a hybrid of ordinary percolation and explosive percolation. However, the potential for tricriticality in other 
correlated percolation models is intriguing and remains a possibility. For example, tricriticality has been explored in 
hybrids of ordinary percolation and explosive percolation in finite-dimensions [28[ . 

Here, we explore the possibility of tricriticality in (1) a mixture of fc = 2 and fc = 3-core edges on random graphs 
and (2) a mixture of fc = 3-core edges and jamming percolation edges in two-dimensions. Cellai and collaborators 
have studied the mixture of fc = 2 and k = 3-core on locally tree-like graphs and random graphs [291 . T herefore, we 
expect agreement between our results obtained using a dynamic rate equation method and theirs [29J. As for the 
second case, jamming percolation models are newer so that we will address some of their properties and argue why 
the search for a tricritical point in such hybrid models may reach a dead end. 

While explosive percolation potentially applies to community networks [3(| and human protein homology net- 
works |3lj |. fc-core percolation may apply to glassy and jamming systems. The Fredrickson- Andersen model |32j is a 
kinetically-constrained model mimicking the caging effect |33[ in glassy dynamics. The onset of a spanning cluster in 
fc-core percolation corresponds to the onset of a glass transition in the Fredrickson- Anderson model [34[ . Recently, the 
Fredrickson-Andersen model was extended to incorporate inhomogeneity in the caging dynamics, i.e. some particles 
are able to break out of their cage if there are less than fc = 2 particles surrounding it, while others require fc = 3 35]. 
This system exhibits a tricritical point in mean-field as the average value of fc is tuned from two to three, the properties 
of which were investigated by Cellai and collaborators [2i| . As for another application of fc-core percolation, the fc-core 
condition encodes the scalar aspect of the principle of local mechanical stability present in jammed packing. It turns 
out that the exponents associated with the mean-field fc > 3-core transition are the same exponents as measured in 
the jamming transition 19]. 

The paper is organized as follows. In the next section (Section II), we present our results for a mixture of fc = 2 
and fc = 3-core on the random graph. In Section III we address a mixture of fc = 3-core and jamming percolation 
models in two-dimensions. In Section VI we discuss the implications of our results. 



II. HYBRID fc-CORE ON RANDOM GRAPHS 
A. Revisiting the rate equation approach 



In 1996, Pittel, Spencer, and Wormald 36] proved that the fc > 3-core transition on random graphs is discontinuous 
using rate equations for the fc-core culling procedure 13711 . This method allows one to easily obtain the discontinuity 
and we shall review it here as has been done in Ref. [38[. 

Consider a random graph with N vertices. At each time step, one vertex whose edges is less than fc is removed 
from the graph. The notion of time is given by time step, t = T/N, where T is the number of algorithmic steps taken 
so far. The change in time is given by At = 1/N, which becomes continuous when N — > oo. At time t, the number 
of vertices is N(T) = (1 — t)N, and the corresponding distribution of connectivity is P z {t) = N Z (T)/N(T). 

One can write down the expected change for N z in the (T + l)th step. It contains two different contributions: 1) 
The first contribution corresponds to the removed vertex itself. It appears only in the equations for N z with z < fc. 2) 
The second contribution is from the neighbors of the removed vertex. The number of its neighbors with connectivity 
z = k will be decreased by 1, hence the total number of vertices with connectivity fc is decreased by 1. The number 
of its neighbors with connectivity z = fc + 1 will be decreased by 1, too, leading to the total number of vertices with 
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connectivity k increased by 1. Therefore, the equation for the change for N z for all z is: 

Here, Xz equals one if z < k and zero otherwise. In addition, \ = XzPz(t) and ~z% — J2 Z z XzPz(t)- The average 
connectivity is defined as Z(t) = z ~J2z z Pz(t)- 

In the thermodynamic limit, the difference equations become differential equations, or 

|((i- t )p.(«)}-^ + fl-^ + «i±^]. w 

This inhnite set of differential equation is not easy to solve. However, one can assume that as vertices get removed, 
those who have never been touched, i.e. whose connectivity is greater than k, obey Poisson statistics with an effective 
connectivity j3(t) and initial connectivity (3(0) = c. In other words, 

{ l-t)P z{t ) = ^P = e-^^f- Vz>k. (3) 
With this ansatz, one can write the normalization condition as: 
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W ) = l-^e-^)^f. (5) 
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Moreover, the average connectivity, Z(t), can be written as: 

oo 



Finally, the rate equation becomes, 



One can also obtain 



z=0 

= ±zXzPz(t) + ±^e-^Mf 

z=0 z=k 

= zx+M-F^it))- (6) 
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(l_ W) = ffl (8 ) 



by comparing $(t) and ^((1 — t)Z(t)) with the latter obtained using the rate equation. 

Now, the culling procedure ends when the graph reaches a stable fc-core configuration. In other words, when 
X = z~X = for all z>k. Therefore, Eq. 5 and Eq. 6 simplify to 

l-tf = F k ((3 f ) 

^ = F k ^(B f ), (9) 
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where /3f = fi(tf) and we have used Eq. 8. To study the nature of the transition, for a given c (or initial occupation 
probability c = p/N), one solves for (3f using the second equation and then computes 1 — tf, which yields the fraction 
of vertices left and is of the order of the giant component should it exist. For k = 2, the critical initial concentration 
signalling the onset of the giant component is c q = 1. For c = c g + e with e « 1, 1 — tf = 2e 2 + 0(e 3 ), i.e. the 
transition is continuous. For k = 3, 1 — tf is finite at the transition such that the transition is discontinuous. 



B. Hybrid model of k = 2-core and k — 3-core 



To search for a tricritical point, we define a model where some fraction of k = k\ vertices, /, and the remaining 
fraction as k — k 2 vertices. In this model, normalization demands that 



ki — 1 oo &2 — 1 

= /( E p ^ + E p ^)) + (! - /)( E p ^ + E p * (*)) ( 10 ) 

z=0 z=fci Z=0 Z=/C2 

As before, we will assume that for the fci-core vertices, as long as z > k\, their Poissonian structure is retained with 
some effective connectivity that changes as a function of time. We will assume this property for the fe-core vertices 
as well with both types of vertices having the same effective connectivity since they are part of the same graph. 
Therefore, 



i = /(X>^) + E i 1 ^^) + (1 - /)(f>^w + E rh^^n)' (") 

z=0 z=fei z=0 z=k2 



where Xfei applies to the fci-core vertices with \k x — for z > k\ and is unity otherwise, and similarly, for Xk 2 with 
&2 replacing k\, to arrive at 

i = fxki+Y^FkAm) + (i - /)0fe + j^FMt)))- (12) 

Moreover, the average connectivity can be written as: 

Z(t) = f{zxkl+ Y$- t F kl -i(J3(t))) + (1 - /)(% + Y^fa-iWt)))- (13) 

We, again, use the rate equation to obtain a relation between the average connectivity and the effective connectivity. 
The rate equation for the hybrid model is 

|«i-»fi<«» - ,^!«[.afc±M w . ffl j^ | .^ + fc±Mj. 

at Xfei Xfci Z(t) Z(t) Xk 2 Xk 2 Z(t) Z(t) 

(14) 

Using the Possionian ansatz, the RHS of the rate equation becomes 



Xk~, (z- 1)1(1 -t)Z(t) z\(l-t)Z(t)' 

+ (1 _ f) MEr_ e -«*) ffi + e -ffW 1 (15) 

1 (z- 1)1(1 -t)Z(t) +e z\(l-t)Z(t) U [Lb) 



and the LHS is 



to arrive at 



«<)[e-«' > g^-e-«' > ^] (W) 



"f^-^w^ 1 -^ (17) 
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where m(t) — (1 — t)Z(t). It turns out that 
dt 



m(t) - l{(l-t)y z P z (t)} 



5>|{(l-t)P,(t)} 



+ /£^[_^ + £(£_!)] _ (1 _ /} £^ + (1 _ /)^[--*L + fcil] 

Xfei Xfci ^(i) Z{t) X k 2 Xk 2 Z(t) Z(t) 



= -2(/= + (1 - /) = ) (18) 

Xfcl X&2 



such that 



/3{t) _ 1 m(i) 



(19) 



whose solution is m(t) = as before. 

Now we are ready to extract the criticial behavior for this hybrid model by occupying edges at random with an 
initial average connectivity Z(0) — c and iterating the culling process until a stable configuration is found at time 
t = tf. Given the fc-core constraints, when t = tf, \k~l = XfcT = z Xk ± = z\k 2 = sucn that the normalization condition 
becomes 

1 = / T -^—e-Pf ?f + y -^—e-Pf P -l (20) 
J i-t f z \ 1-tf z\ y ' 

Z — ki J Z — K2 



with j3f = /3(tf). The average connectivity equation becomes 

k 2 -l 

A 



z(t f ) = f y Y=r e ~ 0f ^\ +pf T=r Fk >- i{M - (21) 



z—ki 



For fci = 2 and = 3, 



Using Eq. 8, we arrive at 



z{tf) = T^V^T + r^7 (1 - e ^ /(1 + M) - (22) 



^ = /e"^ + 1 - e-"' (1 +/?/) = F{fif). (23) 

As before, this self-consistency equation determines /?/ and then one can use the normalization condition at t = tf 
to find the size of the giant component. See Figures 1-3 for a graphical representation of this equation for different 
values of /. 

When / = 1, the model reduces to k = 2-core, and the transition is continuous. When / = 0, the model reduces 
to k — 3-core, and the transition is discontinuous. As / is varied between zero and unity, the curvature of F(f3f) at 
/3f = changes from positive to negative such that for some particular value of /, the curvature of F((3f) at /3f = 
vanishes. In other words, there exists a tricritical point separating the continuous k = 2-core transition from the 
discontinuous k = 3-core transition. This tricritical point occurs at / = 1/2. We will first examine the scaling at the 
tricritical point and then for / < 1/2 and / > 1/2. 

1. The tricritical point: / = 1/2 
When / = 1/2, the self-consistency equation (Eq. 23) reads 

^ = \e-^(3 f + l-e-^(l + f ). (24) 
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FIG. 1: For / < 1/2, the self-consistency equation has two coexisting solutions, /?/ = and finte /3f. The dotted arrow indicates 
increasing c and the bold arrow indicates the transition point 
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FIG. 2: For / > 1/2, the self-consistency equation has one solution for fif. The arrows denote the same as in Figure 1. 



First we examine the scaling of /?/ with c near the transition. The critical value of c indicating the onset of the 
giant component, cg, is determined by F'(j3f)\p f= o — ^ yielding cq = 2 at the tricritical point. Let c = 2 + e with 
e << 1 is a small number and assuming /3f changes continuously (as is indicated graphically) such that /?/ << 1, the 
self-consistency equation becomes 



^(1-- 
2 v 2 



0(e 2 )) 



2 12 yHf) 



(25) 



such that (3f — y/3e^ . 

To find the scaling of the size of the giant component, or (1 — tf)N, as a function of the initial average connectivity, 
c, the normalization condition becomes 
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FIG. 3: For / = 1/2, the self-consistency equation has one solution for /?/, but it scales differently with an increase in the 
initial average connectivity than the / > 1/2 case. The arrows denotes the same as in Figure 1. 



l-U = I e -^ + (l-y e H»/^f 
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Expanding in j3f yields 



i a 1 

l- e -Pt - e -f>tj3 f -± e -PfZL. (26) 



1 - t f = 1 - (1 - /3, + \% - \§ + 0(/3j))(l +/3 f + \$) (27) 



such that 1 — tf = j/3f leads to 



1 - tf = -re. (28) 



This scaling relation is consistent with Ref. |29| |. 

One can also vary / and find the scaling of 1 — tf with /. From Eq. 24, instead of changing c, we change / from 
1/2 to 1/2 + e with < e << 1. After expanding in e and f3f, /3f cx e 1 ' 2 . Since the scaling of (if as both c and / are 
increased beyond the transition, the size of the giant component increases in the same way beyond the transition. 



2. f<l/2 

When / < 1/2, Figure 1 gives us an indication of the scaling behavior near the transition. By increasing the slope 
of the straight line, i.e. decreasing c, the first solution, f3f = (3q > 0, appears when the line is tangent to F(j3f) such 
that the critical average connectivity is 

— = e- p0 [(l-m + f}. (29) 

CG 

Increasing cg by e and assuming/3/ = /3q + Ae x , where A and A are positive constants, and then the LHS of Eq. 23 
becomes 

= — (A) + Ae*)(l - ±) = - + Ae x - (30) 

CG + e CG CG CG CG CG 
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FIG. 4: Increasing / beyond the tricritical value of 1/2, one observes that /?/ changes smoothly. Here, the dotted arrow denotes 
an increase in /. 



and the RHS becomes 



-fh-A* [1 + (1 _ /)(A)+Ae AV| 

1 



- Ae x + -A 2 e 2X - -A 3 e 3X + 0(e 3A ))[l + (1 - /)(#, + Ae x )}. 



6 



(31) 



According to Eq. 30 and Eq. 31, the constant terms and the e x terms on both sides cancel each other, thus, the 
lowest order of the RHS, i.e. e 2A , should cancel the e term on the LHS leading to A = 1/2. To find the scaling of 
the size of the giant component with the initial average connectivity, 1 — tf oc (3f oc B + Ce 1 / 2 , where B and C arc 
positive constants. In other words, the transition is discontinuous. 



3. />l/2 



Figure 2 qualitatively demonstrates that the onset of a nonzero /?/ as a function of c is a continuous transition 
starting at /3/ = 0. The critical value of c is given by 



/--!-. 

CG 

As c is increased by e to cq + e and assuming f3f changes continuously from zero, j3t ~ e. Therefore, 

1 - t f ~ e 2 



(32) 



(33) 



near the transition. This scaling is to be contrasted with the scaling at the tricritical point where the order parameter 
exponent is unity. 

One can also investigate the scaling of the transition with /. As / changes from fo to /o + e, j3f changes from to 
At x . Then, the LHS of Eq. 23 becomes 



Ae x 

c 



(34) 



and the RHS becomes 



1 - (1 - Ae x + ^A 2 e 2X - ^A 3 e 3X + 0(e 3X )) + (/„ - 1 + e)(l - Ae x + ^A 2 e 2X - ^A 3 e 3X + 0{e x ))Ae x (35) 
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FIG. 5: The c, / phase diagram. To the left of the boundary there is no giant component and to the right there is. 



with 

constant term 
e A term 

e 2A term 

e A+1 term 

e 3A term 



The e A terms cancel each other on both sides and e A+1 and e 2A terms sum to zero leading to A = 1 and A — f l _i - 

JO 2 

Thus, the scaling of 1 — tj with a small change in / is the same as for a small change in c. We see that the amplitude 
diverges when /o = 1/2 indicating the vanishing of the /3J contribution such that the cubic contribution comes into 
play when / = 1/2. 

4- Summary 

Our results for the hybrid k\ — 2 and ft2 = 3-core on the random graph can be summarized in the phase diagram 
depicted in Figure 5. For / > 1/2, the k\ — 2-core dominates and the transition is continuous with the size of the 
giant component scaling quadratically with a small increase in the initial average connectivity beyond its critical value. 
For / = 1/2 there exists a tricritical point with a new order parameter exponent, and for / < 1/2 the transition is 
discontinuous. 



III. HYBRID fc-CORE AND JAMMING PERCOLATION MODELS IN TWO DIMENSIONS 

In order to understand what happens in the hybrid models, we first review two jamming percolation models since 
they are rather new to the field. These models are vertex models (no edges are explicitly added). Of course, a pure 
vertex version of fc-core percolation can also be introduced where an occupied vertex requires at least fc occupied 
vertices to remain occupied. It is this version of fc-core we will refer to in the following section. 
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A. Counter-balance model 

While k > 3-core exhibits a discontinuous transition in mean-field, there is no known two-dimensional example of 
k > 3-core exhibiting a discontinuous transition. In fact, it appears that k > 3-core on two-dimensional lattices either 
exhibits an ordinary percolation transition with a shift in the critical occupation probability, or there is no transition 
until the lattice is fully occupied [2(| Hl| . For instance, simulations of k — 3-core percolation on the triangular lattice 
result in similar ordinary percolation exponents, while for 4 > k > 6, p c = 1. A heuristic argument behind the former 
result is that finite clusters are allowed for k = 3 with a fully occupied hexagon being the smallest structure. One can 
then imagine this object to be smallest object, as opposed to a single vertex, such that path-like spanning structures 
are formed out of fully occupied hexagons. This procedure is a simple coarse-graining on a microscopic scale and will 
not effect the macroscopic scales near a continuous phase transition. 

Inspired by the jamming transition in two-dimensions where the fraction of particles participating in the jammed 
structure goes from zero to finite [39l |40T| . one can encode various properties of the jammed packings into a percolation 
model. One of those properties is counter-balancing. For the force on each particle to be balanced, there must be 
an occupied particle on either side. Of course, in two-dimensions, two particles on opposite sides of the particle in 
question will not suffice since the configuration is not mechanically stable. However, three particles whose centers are 
120 degrees with respect to each other as measured from the center particle is a stable configuration. 

Counter-balance percolation takes into account the counter-balancing aspect of force-balance [4l| . As for an exam- 
ple, we begin with a two-dimensional square lattice. Each vertex neighbors all vertices within a 5x5 square modulo 
itself. In other words, each vertex has 24 nearest neighbors. The counter-balancing constraint is the following: for 
an initially occupied vertex to remain occupied, there must be at least one occupied neighbor in set A, which in turn 
calls for at least one occupied neighbor in set B, and there must be at least one occupied neighbor in set C, which 
in turn calls for at least one occupied neighbor in set D. The four sets A, B, C, and D, are defined in Fig. 6. The 
counter-balance constraint can be succinctly stated as: (A and B) and (C and D), where each letter X is short for 
"at least one occupied vertex in set X" . Note that the counter-balance constraint is defined in such a way such that 
vertical and/or horizontal lines of occupied vertices are, by themselves, not stable. Fig. 7a demonstrates an allowed 
configuration and Fig. 7b demonstrates a forbidden configuration. 

To enforce the counter-balance constraint, we initially occupy vertices on the lattice with independent occupation 
probabilities p, and then repeatedly remove occupied vertices that violate the counter-balance constraint, until all 
remaining occupied vertices obey the constraint. Note that p is the occupation density before culling, and generically 
differs from the final occupation density. Moreover, the model is abelian, i.e. the order of the culling does not affect 
the final configuration. 





















A 




















B 
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FIG. 6: Force-balance model on the 2d square lattice with 24 nearest neighbors. 



Numerical simulations of the counter-balance model strongly suggest a discontinuous transition with a finite fraction 
of vertices participating in the spanning cluster at the transition [4lj . Of course, numerical simulations can be 
misleading as has been demonstrated in the past, and so one looks for evidence beyond numerics. Such evidence is 
provided for by studying a simpler, but related, jamming percolation model, namely the spiral model. While the 
spiral model is less physical, one can make several concrete statements about its percolation transition. 



11 



(a) 



(b) 



FIG. 7: Black sites denote occupied sites, and red sites denote unoccupied sites, (a) Allowed configuration, (b) Forbidden 
configuration. While the number of occupied nearest neighbors is greater than three, the force-balance condition is violated. 



B. Spiral model 

The spiral model [I?} is defined as the following. Again, we begin with a square lattice. The neighbors of each 
vertex contain the four nearest neighbors and the four next nearest neighbors to make a total of eight neighbors. After 
the initial random occupation of vertices, for each vertex to remain occupied there must be at least one occupied 
neighbor in set A and at least one occupied neighbor in B, or there must be at least one occupied neighbor in set C 
and one occupied neighbor in set D. See Fig. 8 to denote the sets. 

It turns out that the critical occupation probability for this model is the same as for directed percolation [421 ] . 
denoted as p® p . One can see the link with directed percolation when considering only sets A and B or sets C and D. 
Each pair of sets is isomorphic to the canonical two-dimensional version of directed percolation with two neighbors 
above and below each vertex. Therefore, for p > p® p , there exists a spanning cluster along either diagonal. Since 
there are two pairs of sets one might argue that the critical occupation probability is lower than that of directed 
percolation. However, one can show that a certain class of voids (finite "clusters" of unoccupied vertices) cause 
all remaining occupied vertices in the system to become unoccupied as long as p < p® p . Given the two bounds, 
Pc P — P? P i where SP denotes the spiral model. 
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FIG. 8: The sets in the spiral model. 



While determining p~ p is a detail, it is an important one for arguing that the percolation transition is discontinuous, 
which is very different from the continuous transition in directed percolation. To argue for a discontinuity, one can 
construct a set of spanning structures, which include the origin, and demonstrate that for p = p® p , the probability of 
such a set is greater than zero. In other words, the spanning structure is compact at the transition. This argument 
involves the notion of T-junctions, which characterize one diagonal path being supported by the other diagonal path. 
For instance, a nonspanning path along one diagonal can survive only if it is sandwiched between two paths of the 
other diagonal. See Figure 9. 

To observe compact spanning structures at the transition, one can build a spanning structure beginning with a 
rectangle whose long side is along the A-B diagonal and includes the origin. Then, one considers an infinite sequence 
of pairs of rectangles of increasing size emanating outward from the initial rectangle and intersecting as indicated in 
Figure 10. If each of the A-B rectangles contains an A- B spanning path along its length, it is indeed a spanning 
structure containing the origin due to the existence of T-junctions. From what is known about directed percolation 
at the transition, one can show that the probability for each of the A-B rectangles to contain a A-B spanning path is 
indeed greater than zero such that this compact spanning structure exists at the transition. Please see Ref. [27] for 
details. 

It turns out that one can extend these arguments to jamming percolation models with more than two pairs of 
sets (43j . However, these arguments have not yet been extended to sets with more than two sites. The universality 
of directed percolation should allow for such an extension since as long as the sets are arranged in opposite pairs, 
the occupation for each pair of sets is governed by a directed percolation-type process. It turns out that the rules 
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FIG. 9: A pair of T-junctions in the spiral model. 



for the counter-balance model can be written as a spiral-like model with more than two pairs of sets with some sets 
containing more than two vertices. See Figure 11. In addition, the conversion also involves triplets and quadruplets of 
sets (as opposed to just pairs) with each set participates in, for example, a pair-wise interaction as well as a three-way 
interaction. These triplets and quadruplets of sets have not yet been addressed in the context of a spiral-type model 
but should not invalidate the overall construction of the above argument. The fact that one set participates in several 
interactions increases the possible configurations, but, again, should not invalidate the above argument. Therefore, we 
expect the percolation transition in the counter-balance model to be discontinuous with numerical evidence supporting 
this expectation [411 ] . 



FIG. 10: Schematic of a spanning structure in the spiral model where a path spans each rectangle. 



D' 



A' 



H' 



FIG. 11: The counter-balance rules can be implemented with (A' and C) or (D' and B') or (E' and B' and C) or (F' and C 
and D') or (G' and A' and B') or (H' and A' and D') or (A' and G' and H') or (B' and E' and G') or (C and E' and F') or 
(D' and F' and H') or (E' and F' and G' and H'). 
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C. Tricriticality in the hybrid fc-core/spiral model? 

For the spiral model, the notion of a T-junction as well as properties of directed percolation are key in establishing 
a discontinuous transition. What happens when the model is perturbed by having an infinitesimal fraction of vertices 
with the fc-core condition replacing the spiral condition? To answer this question, we must first choose fc. The value of 
fc is chosen such that p c for the lattice with every vertex obeying the fc-core condition is less than p® p . For example, 
k = 3 satisfies this condition. Note that each vertex has the same eight neighbors as in the spiral model. In addition, 
when all vertices are fc-core vertices, we assume that the transition is continuous and in the same universality class as 
ordinary percolation. This model is similar to the fc = 3-core transition on the triangular lattice. In both cases, there 
exist small clusters that survive the culling process such that, heuristically, these clusters link up to form path-like 
structures as occurs in ordinary percolation. 

Given the properties of the all fc-core transition, as the fraction of spiral model vertices increases from zero, the 
critical occupation probability increases since the fc-core condition is less constraining than the spiral model condition. 
As long as the critical occupation probability is less than p^ p , the construction invoking a spanning scaffold of T- 
junctions at p® p to demonstrate a discontinuous transition for the all spiral model vertices no longer holds. The 
integrity of the T-junctions to arrive a spanning structure is destroyed with finite clusters now allowed even for an 
infinitesimal fraction of fc-core vertices. T-junctions are no longer necessary to support a structure. This property is 
consistent with the violation of the former bound, p > p® p ' , since each diagonal is no longer isomorphic to directed 
percolation independently. The second former bound, p < p® p , also breaks down since the growth of voids is now 
stopped by fc-core vertices. 

The existence of finite clusters for any fraction of fc-core vertices is certainly different from the spiral model where 
no finite clusters are allowed. Does the existence of finite clusters imply that the transition is continuous up until 
all fc-core vertices are replaced with spiral model ones? In fact, there's no direct relationship between the existence 
of finite clusters and the continuity of a transition. One can construct a model occupying vertices at random and 
independently and then remove all finite clusters. While this particular constraint is highly non-local, it preserves the 
properties of the spanning cluster (of ordinary percolation) at the transition. Even so, we argue that the transition 
is continuous for the following reason. As the fraction of spiral model vertices increases, their increasing presence 
demands the increasing use of T-junctions to support spiral model vertices where two paths along one diagonal 
sandwich and support the path along the second diagonal. However, these two supporting paths can now end on 
fc-core sites as opposed to other T-junctions to form finite clusters. Two finite clusters whose interiors each contain 
spiral model sites can join in such a way that the removal of one site in the two cluster formation that is not shared by 
both clusters before the joining does not induce the removal of the other cluster. We speculate that this independent 
cluster joining property leads to a continuous transition since the joining of such clusters leads to path-like structures 
on a larger scale as clusters are placed "side-by-side". The larger the clusters, the larger the scale one has to go to 
"observe" the path-like structures. See Figure 12. 

One would like to make the speculation that the independent cluster joining property leads to a continuous transition 
more rigorous. Such speculation may lead to a framework to prove that the transition for similar models, such as 
fc = 3-core on the triangular lattice, is in the same universality class as ordinary percolation. Currently, there is only 
numerical evidence for fc = 3-corc on the triangular lattice being in the same universality class as ordinary percolation. 

o o o o o o o 
o o • o • o o 



o ••[•]•• o 

o o • o • o o 

o o o o o o o 

FIG. 12: An example of two finite clusters for k — 3-core in the eight neighbor model sharing a boxed occupied vertex. The 
removal of any occupied vertex to the right of the boxed occupied vertex does not trigger the removal of any occupied vertex 
to the left of the boxed occupied vertex. 



So, even for an infinitesimal fraction of fc-core vertices with all other vertices dictated by the spiral model, we 
speculate that the transition is in the same universality class as ordinary percolation. In other words, only when all 
vertices are spiral model vertices is there a discontinuous transition. Given this scenario, we speculate that there is 
no tricritical point in this hybrid of fc = 3-core and the spiral model. 
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FIG. 13: Differential for the probability of spanning, P S p an , as a function of occupation probability p for g = 0. Different 
system sizes are studied. The inset demonstrates the finite-size scaling for P spa n(p) using the correlation length exponent for 
ordinary percolation. 



D. Tricriticality in hybrid fc-core/counter-balance model? 



Since the counter-balance model can be expressed as spiral model-type constraints, what happens in the hybrid 
fc-core/spiral model may resemble what happens in the hybrid fc-core/counter-balance model. Therefore, we expect 
that the transition is continuous as long as 1 — / = g < 1, where g is the fraction of counter-balance vertices, with 
fc = 3. We now provide numerical evidence to support this claim. 

First, in Figure 13 we plot the differential curve for the probability of spanning, P sp am as a function of p for 
g = and fc = 3 with 24 nearest neighbors on the square lattice. The system length is denoted by L and periodic 
boundary conditions are implemented. The position of the peak denotes the critical occupation probability, p c , for 
each particular system size. We then perform a scaling collapse using the correlation length exponent of ordinary 
percolation v = 4/3. See the inset to Figure 13. The scaling collapse is reasonable suggesting that the transition is at 
least consistent with the ordinary percolation universality class as should be the case for g < 1 given the discussion 
in the previous subsection. 

Next, we plot average size of any spanning cluster, < S > (the brackets denote the configuration averaging) as a 
function of the initial occupation probability p for different gs. See Fig. 14. In an infinite system, for g < 1, we expect 
< S >— A(g)(p — p c (<7)) 5 / 36 for p > p c {g) and zero otherwise. In a finite system, this function will be smoothed out. 
As expected, we observe that p c (g) increases with increasing g. The curves, even for g = 1, appear to be sharpening 
with increasing g as well such that the amplitude A(g) increases with increasing g. However, the < S > versus p 
curves appear qualitatively similar, so how does one differeniate between a continuous transition and a discontinuous 
one? 

We can do so by examining the distribution of the sizes of the spanning cluster as opposed to just looking at the 
average size. See Figs. 15 and 16. To compare distributions among the different gs, we plot the distribution for the 
same average spanning size, which means p changes from curve to curve. The distribution for g — has a well-defined 
peak at small values of S, small meaning much less than unity. As g increases from zero, the distribution for S looks 
similar with a slight shift in the distribution. It is not until g « 0.99 that the position of peak increases such that it 
is much greater than the fixed average size of 0.05 indicating that in many instances no spanning cluster found and 
when they are found, the spanning cluster is much larger than the average. This feature of the distribution, as long 
as it persists in the infinite system limit, is characteristic of a discontinuous percolation transition. When the system 
size is increased for fixed g, we observe that the position of the peak shifts to the left indicating that the transition 
becomes more continuous even for g = 0.99. Therefore, the numerical data supports our scenario of the absence of a 
discontinuous transition until g = 1 and no tricritical point. 
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FIG. 14: Average size of spanning cluster, < S >, as a function of p for a mixture of k = 3-core and counter-balance percolation 
vertices. For g = 0, the vertices are all fc-core. For g = 1, the vertices are all counter-balance. Even for g = 0.96, the curve 
resembles the g — (continuous) case. Here, the system length, L, is L = 256. 



IV. DISCUSSION 

There has been a recent proliferation of percolation models going beyond the original model, where edges/vertices are 
randomly and independently added to a network, to find discontinuous transitions. While much of this proliferation 
is driven by what is called explosive percolation, there has been, shall we say, quieter progress on what is called 
jamming percolation models inspired by jamming and glassy systems. Interestingly, there exists at least one example 
of a discontinuous percolation transition on random graphs, called fc-core percolation, which was posed back in 
1979 [ItI [36j . However, this transition is not the ordinary type of discontinuous transition lacking any diverging 
lengthscales. In fact, the k > 3-core transition exhibits several diverging lengthscales. In addition, the spiral model, 
and its related counter-balance model, provide examples of discontinuous percolation transitions in two-dimensions, 
again, with diverging lengthscales [13, |4l[ . 

Since both k > 3-core and the spiral model exhibit discontinuous transitions, it may be useful to investigate what 
properties they share in order to execute a more efficient search for other discontinuous percolation transitions. Both 



0.07 




FIG. 15: Distribution of the size of spanning cluster, S, for different gs. For g = 1, the vertices are all counter-balance. Only 
for g = 0.99 does the distribution begins to shift to the right significantly, whereas for g — 0.98 the distribution more closely 
resembles the g = case. Here, again, L — 256. The no spanning cluster contributions are not shown. 
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FIG. 16: Distribution of the size of spanning cluster, S, for different gs and for different system lengths. Note that as L 
increases, the curves for g — 0.99 shift to the left indicating that as L becomes large, the transition approaches the g — 
continuous case as opposed to the g = 1 discontinuous case. 



models share the property of not "allowing" finite clusters. In the spiral model, this property is exact. For the 
fc-core model on trees, this property is exact. On random graphs, this property is approximate but becomes an 
increasingly better approximation in the thermodynamic limit. One may argue that forbidding finite clusters implies 
a discontinuous transition, at least in low-dimensions. However, this is not the case, since one can construct the usual 
percolation model on the square lattice and remove any finite clusters using a non-local rule. One can also do this 
using a local rule using only one pair of sets from the spiral model. 

The nontrivial nesting of the spanning structure in the spiral model with the two diagonals interdependent on one 
another tells us that even if finite clusters are allowed, the formation and subsequent joining of even finite clusters 
should be non-local in the sense the removal of one vertex in a cluster signals the removal of at least a finite fraction 
of vertices in its formerly disjoint neighboring clusters. The rarefication [44j of the lattice or a cluster aggregation 
model with a non-local kernel [45|, |46| provide such mechanisms for nontrivial cluster joining such that the transition 
is discontinuous. On the flip side, we speculate that perhaps the independent finite cluster joining property, where 
the removal of one vertex in one of the two finite clusters (determined prior to joining) that is not shared by the other 
does not induce the removal of the other finite cluster. This property provides for a simple path-like joining of finite 
clusters such that an ordinary percolation transition presumably occurs as is found in k > 3-core models in finite 
dimensions. 

Given the existence of both continuous and discontinuous percolation transitions, one may form hybrid models of 
the two types to search for a tricritical point separating the continuous regime from the discontinuous regime. Using 
a mixture of k — 2-core and k = 3-core vertices on random graphs, we find a tricritical point and determine the 
size of the giant component as a function of the average connectivity. While two previous works have identified this 
point [29l . [43] , we use a different dynamical method via a rate equation approach, which could prove to be valuable 
for investigating other correlated percolation models as the constraints become complex. Of course, our results using 
the rate equation approach agree with previous results using a "static" method. 

Moreover, we investigate the possibility of a tricritical point in two-dimensions. We do so with a mixture of 
k = 3-core vertices and counter-balance vertices where the full k = 3-core model exhibits a continuous transition 
(which differs from mean-field) and the full counter-balance model exhibits a discontinuous phase transition. We 
argue that there is no tricritical point in this mixed model since the discontinuous transition occurs only in the full 
counter-balance model. However, we expect interesting crossover behavior between to the two types of transitions that 
should be explored. This result is the first we know of with a line of continuous transitions ending at a discontinuous 
transition, which is to be compared with the water phase diagram where there is a line of discontinuous transitions 
ending at a continuous one. In addition, this result is to be contrasted with the work of Cellai and collaborators 
studying a mixture of k — 2 and k = 3-core vertices on the square lattice 29]. For the full k = 2-core model on 
the square lattice, the transition is continuous, while for the full k = 3-core model, p c = 1 with interesting crossover 
behavior between the two cases. This result is also to be contrasted with the tricriticality obtained in a diluted Q-state 
Potts model on a triangular lattice [48| . 

Finally, what about ties between mixed correlated percolation models and physical systems? The introduction 
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of hetereogeneities into the Frederickson- Anderson model for glassy systems maps to the mixed fc-core model [35| . 
Moreover, there is a recent experiment [4^ | where a two-dimensional packing of frictional disks is sheared at a packing 
fraction just below the shear-free jamming transition to induce jamming. At small applied shear stress, a subset 
of the contact network of the jammed states exhibits spanning structures along one direction only, while at larger 
applied shear stress, the spanning structures (of a subset of the contact network) percolate in both directions. One 
can perhaps model the connectivity of the spanning structures as the applied shear stress is varied by taking a mixture 
of the usual counter-balance model and an anisotropic version of the counter-balance model where only sets A and 
B are considered. More specifically, as the applied shear stress is increased the ratio of counter-balance vertices to 
anisotropic counter-balance vertices increases. However, since only a subset of the contact network is used in obtaining 
this result, some finite structures may be allowed. 

JMS acknowledges support from NSF-DMR-CAREER Award 0645373. 
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